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ABSTRACT 


We amplify previous discussions of the fine-tuning price to be paid by supersymmetric models 
in the light of LEP data, especially the lower bound on the Higgs boson mass, studying in par¬ 
ticular its power of discrimination between different parameter regions and different theoretical 
assumptions. The analysis is performed using the full one-loop effective potential. The whole 
range of tan/3 is discussed, including large values. In the minimal supergravity model with 
universal gaugino and scalar masses, a small fine-tuning price is possible only for intermediate 
values of tan/3. However, the fine-tuning price in this region is significantly higher if we require 
b — t Yukawa-coupling unification. On the other hand, price reductions are obtained if some 
theoretical relation between MSSM parameters is assumed, in particular between /x 0 , Mi/ 2 and 
Aq. Significant price reductions are obtained for large tan/3 if non-universal soft Higgs mass 
parameters are allowed. Nevertheless, in all these cases, the requirement of small fine tuning 
remains an important constraint on the superpartner spectrum. We also study input relations 
between MSSM parameters suggested in some interpretations of string theory: the price may 
depend significantly on these inputs, potentially providing guidance for building string models. 
However, in the available models the fine-tuning price may not be reduced significantly. 
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1 Introduction 


To paraphrase Saint Augustine: ’’May Nature reveal supersymmetry, but not yet.” This seems 
to be the message from LEP and other accelerator experiments, so far. There are tantaliz¬ 
ing pieces of circumstantial evidence for supersymmetry at accessible energies, including the 
measured magnitudes of the gauge coupling strengths [||, |J and the increasing indication that 
the Higgs boson may be relatively light |3|, f|]. On the other hand, direct searches at LEP 
and elsewhere have so far come up empty-handed. In the case of LEP 2, the physics reach for 
many sparticles has almost been saturated, though the direct Higgs search still has excellent 
prospects. 

In this context, it is natural to wonder whether the continuing absence of sparticles should 
disconcert advocates of the Minimal Supersymetric Extension of the Standard Model (MSSM). 
After all, the only theoretical motivation for the appearance of sparticles at accessible energies 
is in order to alleviate the fine tuning required to maintain the clectroweak hierarchy ||, and 
sparticles become less effective in this task the heavier their masses. Since the problem of 
fine-tuning is a subjective one, it is not possible to provide a concise mathematical criterion 
for deciding whether enough is enough, already. Moreover, the fine tuning can be discussed 
only in concrete models for the soft supersymmetry breaking terms, and any conclusion refers 
to the particular model under consideration. The fine-tuning price may also depend on other, 
optional, theoretical assumptions. 

The idea which we prefer to promote here is that, along with the overall increase in the 
fine-tuning price imposed by the data, any sensible objective measure of the amount of fine 
tuning becomes an interesting criterion for at least comparing the relative naturalness of various 
theoretical models and constraints, and - within a given framework - of different parameter 
regions. 


We have recently shown that the latest LEP and other data which constrain the MSSM 
parameters significantly increased the requisite amount of fine tuning [||] compared with pre- 
LEP days. We used one particular measure ]7|, |[, namely A 0 = max^/Mf (<9Mf/<9aj)|, 
where the a% are the input parameters of the MSSM (for other measures of fine tuning, see 
@). Our tree-level analysis clearly demonstrated several qualitative trends but, as an obvious 
improvement, one should use the best available theoretical tools to evaluate the fine tuning, 
including in particular the full one-loop effective potential of the MSSM JTT], |TT|. Secondly, 


one should update the analysis with the most recent experimental information, in particular 
on the mass of the Higgs boson El and the new result for BR(b —> sy) 


With the above improvements in hand, in this paper we address anew the question of the 
necessary amount of fine tuning, with a particular view to the power of the fine-tuning price to 
discriminate between different parameter regions and different theoretical assumptions. 


In Section 2 we recall our measure of fine tuning and discuss its various qualitative aspects 
in the supergravity-nrediated scenario with universal gaugino and scalar masses (the minimal 
supergravity model). The particular role of the Higgs boson mass is elucidated. In Section 3 
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we present our full one-loop results for small and intermediate tan /? in this model. In the first 
place, we confirm previous findings |L0], [11] that including the full one-loop effective potential 
reduces the apparent amount of fine-tuning by about 30% at moderate tan/? ~ 10, and by 
much larger factors for both small and large tan (3. On the other hand, the latest experimental 
lower limit Mh > 90 GeV for low tan /3 increases the price again, so that the fine-tuning price 
we find for low tan (3 is not very different from that in |6j. In the minimal supergravity model, 
for intermediate tan (3: 3 < tan/3 < 15, there still exist domains of the parameter space with 
moderate, 0(10%), fine-tuning. This result is obtained after including all available experimental 
constraints, including in particular b —> sy, but with no constraint on the Yukawa sector. 

In Section 4 we address the question of bottom-quark/tau Yukawa-coupling (b — r) unifica¬ 
tion |L4] in the minimal model for small and intermediate tan (3 ^ 30. We show that inclusion 
of one-loop corrections to the bottom-quark mass substantially enlarges the tan/? region where 
b — t unification is possible, albeit at the expense of a higher fine-tuning price. Furthermore, 
the interplay of the constraints from b — r unification and b — ► sy decay and the dependence 
on tan (3 is understood. The minimal model with both b — r unification and the b —> sy con¬ 
straint imposed has no regions of small fine tuning. However, it is stressed that b —> sy decay 
is an optional constraint, which can be relaxed if we admit some departure from the minimal 
model, e.g., some flavour structure in the up-squark mass matrices. Given such a generalization, 
regions with low fine tuning exist with or without b — r unification. 


In Section 5 we emphasize the dependence of the fine-tuning price on the choice of the 
set of independent soft mass parameters in a given model. In particular, we find that in the 
minimal supergravity model A 0 may be significantly reduced if the parameters /i 0 and M 1 / 2 
or A 0 (depending on the value of tan/?) are considered as linearly dependent on each other. 
In some stringy models these parameters are indeed not independent, although the correlation 
may not be linear. As is briefly discussed in Section 6, we find that in one class of such 
models Ao may be minimized only in unphysical regions of the parameter space corresponding 
to small sparticlc masses and/or the absence of electroweak symmetry breaking. Within the 
physical region of parameters in the models studied, the fine-tuning price may not be reduced 
significantly. 


In Section 7 we discuss the case of large tan /? > 30. Our main conclusion is that it remains 
attractive (with small fine-tuning) for non-universal Higgs boson mass parameters at the GUT 
scale. Section 8 contains our conclusions. 


2 Measure of Fine Tuning and Tree-Level Discussion 


We first specify more precisely the fine-tuning criterion we use. Following |7|. Q, we consider 
the logarithmic sensitivities of M z with respect to variations in input parameters ay 


ai dM z 
M z da 
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Note that here we take derivatives of and not, as in [b]. of M§: hence onr A a . are smaller 
by factors 2, other things being equal. We then define 


An = max I A 


O'i I 


( 2 ) 


It is clear that the fine tuning can be discussed only in concrete models for the soft super¬ 
symmetry breaking terms, and with a specified scale for their generation. Calculation of the 
derivatives (jl]) requires minimization of the effective scalar potential written in terms of the 
a*. In the first approximation one may use (as we did in our previous paper) the tree-level 
form of the potential, but it is known [Rj and has been strongly re-emphasized recently (TT 


that reliable quantitative analysis requires use of the full one-loop effective potential. This 
is particularly important for low and large values of tan/3, where one-loop corrections to the 
tree-level potential are decisive for electroweak symmetry breaking. In this paper we follow the 
one-loop approach of [[HJ. 

It is, nevertheless, useful to discuss first certain qualitative features of our analysis starting 
with the tree-level potential: 


V = rri'i | H] | 2 + m 2 \H 2 \ 2 + U/j(e//,// 2 + eH x H 2 ) 
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The derivatives (|l|) then read 
2 


A a . = 


(t 2 _ 1)2 E |(44 + 44)^ 2 + f) 


didj didj 


Ml 


-c'{ (f+1) 


2 -| \ Q*idj 


Ml 


+ 2 1 


2 d^dj 


M\ 


Ml 


- 4 i 2 


7.2 i -) \ a i a j _|_ 2°*^ 


0 t 2 + 1 ) 


Ml 


Ml 


(4) 


where t = tan/3 and the coefficients Ci 2 , and are defined by 

m l,2 ( M z ) = clzdidj c'l 2 = cf 2 

ij 

s ( M i) = E4«i 


/i 


K) = Ek 


(5) 


The numerical values of the coefficients c{ 2l c l B and c) t can be found by solving the renormalization- 
group (RG) equations [|T5| [Rj for the running from some initial scale down to Mz- 


The most popular model, with some phenomenological backing, is the MSSM with universal 
gaugino and scalar masses (Mi/ 2 ,mo) at the input supergravity scale (or GUT scale), and a 
universal trilincar (bilinear) supersymmetry breaking parameter A$(Bq). The model is then 
formulated in terms of these four parameters and the /io parameter. Important features of this 
model are strong correlations between soft terms: all scalar mass parameters are assumed to be 
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equal. At present we do not have any convincing theory of soft terms and might equally well 
contemplate the possibility of other patterns for them, for instance of non-universal Higgs-boson 
masses and/or different sets of independent parameters. The amount of fine tuning depends 
on this choice, as discussed in Section 5. However, in this and the following Section we discuss 
the universal case with five independent parameters (the minimal supergravity model). 


Several qualitative effects in the fine tuning of soft terms can be seen already from (|4|). In 
particular, we can discuss the typical magnitude of the derivatives taken with respect to the five 
parameters of the minimal supergravity model, with the scale of the generation of soft terms 
taken to be Mqut = 2 x 10 16 GeV. As an example, we consider the region of small tan/3, not 
far from the quasi-infrared fixed-point solution for the top-quark Yukawa coupling, in which 
the fine tuning is generically larger than for intermediate values of tan /3. 


Using analytic solutions obtained in |15| for the coefficients c/ 2 , c l B and c/ as an expansion 
in the parameter 


V = Yt/Y t 


FP 


( 6 ) 


where Y t FP is the fixed-point value for the top-quark Yukawa coupling Y t = h 2 / 47r, we can 
calculate the derivatives in eq. (|4|) explicitly. In the limit y 1 one gets: 


A/i0 (f - l) 2 




£ _ 

Ml 


+ t 2 4 


JY_ 

M\ 


- 1 - 


MX 

Ml 


^ W/2 ~ 


t 2 +1 
(t 2 -1) 2 


(jt 2 - 1) 


M 


1/2 | t ( ^Mi/ 2 ^ yMy 


Ml 


V Ml 


+ 


M\ 


+ 


12f 2 M 2 /2 
t 2 + 1 M 2 a 


A 


Ao 


t(t 2 + 1) 
(f 2 — l ) 2 


/^A 0 /iA 0 

Ml + Ml 


(7) 


A _ 2t( l + !) (l lB 0 »Bo\ 

Bo ~ (t 2 - l) 2 V Ml Ml) 


A 

L -*mo 


(t 2 - 1) 2 


(t 2 + i)(t 2 


2 ) 


Ml 


-21 2 


™ o 

Ml 


where we may consider the region t = tan f3 ~ (1.5 — 2) to be consistent with our approximation 
y 1 - qualitatively similar conclusions can be drawn for other values of tan /3, as long as the 
bottom quark Yukawa coupling is much smaller than Y t . The parameter /i at the scale Mz is 
related to its initial value /io by the equation y 2 ~ 2yl(l — y ) 1 / 2 . We note that the largest 
derivatives are A /i0 and Am 1/2 , and they are of opposite signs. For instance, for all parameters 
of order Mz (a situation already strongly excluded by the present experimental constraints) 
both are already greater than ~ 10 for tan/3 ~ 1.5. They increase quadratically with the values 
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of the parameters, and this is the reason why large fine tuning is found for low tan (5 within 
the experimentally-acceptable parameter range. The derivative A^ 0 is also sizeable and may 
also play an important role, since large negative A 0 may be necessary |1| to satisfy the present 
Higgs-boson mass limit [|T2|. The derivative A mo is typically of little importance. For a given 
parameter set, the necessary fine tuning is determined by the maximal derivative Any 
such set should be chosen consistent with the present experimental data and the constraints 


(correlations) imposed by proper electroweak symmetry breaking (see for instance [|T0 


Among the experimental constraints, a special role is played by the Higgs-boson mass limits. 
This effect can be isolated by imposing all the available experimental constraints except the 
lower limit on Mh■ For a chosen value of tan/3, we get then some minimal value of Ao, and the 
corresponding parameter set determines the mass of the Higgs boson 
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where f(x,y ) = 2 — {x + y)/{x — y) log (x/y). We display in (||) only the one-loop formula valid 
for Ma ~ 3 Mz jl7|, which is a good approximation to the two-loop RG-improved prescription 
[ |T§[ | used in our numerical studies. Due to the logarithmic dependence of M( on the physical 
stop masses squared (which in turn are functions of the initial parameters) any departure from 
the “best” value of Mh, for fixed tan/3 and in the range allowed by the other constraints, 
transmits itself into an approximately exponential rise of A 0 . This happens for M ( changing 
in both directions, towards values both smaller and larger than the “best” value. Thus, for a 
given tan /3, the Higgs boson mass is a crucial probe of fine tuning. We also recall [l8j that the 
Higgs boson mass is maximal for large \A t \ = \A t — y cot (5\. Since A t is given by [0 


A, « (1 - y)A 0 - 0(1 - 2)M 1/2 , (9) 

maximizing Mh requires y > 0 and large negative Ao- Hence the derivative A^ may be large 
in the low tan /3 region. 


3 Results for Low and Intermediate tan [3 


In this Section we discuss our full one-loop results in the minimal supergravity model with 
five independent parameters. It is appropriate to begin by re-emphasizing that the fine-tuning 
criterion is not a rigorous mathematical statement, but rather an intuitive physical preference 
and hence remains necessarily subjective. (For instance, as already mentioned, our present 
definition differs by a factor 2 from the one used in [|6|.) Nevertheless, for a chosen measure 
of fine-tuning, one can study in this model relative changes in the amount of fine tuning as 
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a function of changing experimental limits and of the considered parameter range. Here we 
emphasize this use of the naturalness criterion. 


We first recall the experimental constraints used in this analysis. The data we take into 
account include the precision electroweak data reported at the Jerusalem conference [|4| , which 
are dominated by those from LEP 1. We constrain MSSM parameters by requiring that Ay 2 < 4 
in a global MSSM fit [^1], |22], [2/|, ^4|. The main effect of this constraint is a lower bound on 
the left-handed stop: Mf L > 300 — 400 GeV |23|]. We also take into account the direct LEP 2 
lower limits on the masses of sparticles |^5[ and Higgs bosons. For the latter, we base ourselves 
on the recent data reported in |L2|. We use the limit M h > 90 GeV which, strictly speaking, 


is valid for the Standard Model Higgs boson. This is approximately valid also for the MSSM 
Higgs boson for small tan/3, although there still exist small windows in the parameter space 
where the experimental limit is lower. We neglect this possibility in the present analysis. 


The final accelerator contraint we use is the recently-measured value of the b 
ing ratio 2 x 10 -4 < B(B —> Al,y) < 4.5 x 10 -4 


at the 95% C.L. ]T3 


sy branch- 
The interpretation of 

this measurement in the MSSM is still subject to some uncertainty, because not all the 0(a s ) 
corrections have yet been calculated. Resumming large QCD logaritms of the type log (M w /m^) 
up to next-to-leading order (NLO) accuracy has recently been accomplished |26|. These cal¬ 
culations are identical in the SM and the MSSM. The initial numerical values of the Wilson 
coefficients at the scale /r ss M\y are, however, different in the two models. In our analysis we 
have used for them two-loop results available for the standard W^t and H ± t |27|] contributions, 
and only the leading-order results for the chargino-stop contribution |28| . The uncertainty due 
to 0[oL s j 7r) corrections to them has been, however, included as in |29], [23]. Those references 
also contain extensive discussions of the role played by the b —> sy measurement in constraining 
the parameter space of the MSSM. 


An important role may also be played by non-accelerator constraints, in particular the relic 
cosmological density of neutralinos N°, if these are assumed to be the lightest supersymmetric 
particles, and if R parity is absolutely conserved. Both of these assumptions may be disputed, 
and a complete investigation of astrophysical and cosmological constraints is beyond the scope 
of this analysis (for steps in this direction, see pOj). 


One-loop corrections to the effective potential are taken into account as in [10;]. using the 
decoupling method of |3T]. Numerical calculation of the derivatives A ai is also explained in 
[lCH . Electroweak symmetry must, of course, be broken, and this requirement imposes strong 


constraints on the allowed parameter region. The main effect is that /i > max(My 2 , itlq). 


In Figs. 1 to 3 we show Ao = max|A a J as a function of some mass parameters and some 
physical masses, for tan (3 = 1.65, 2.5 and 10. The results are in agreement with the qualitative 
discussion of Section 2 and with the results of 0 , but the inclusion of one-loop corrections to the 
scalar potential sizeably decreases the fine-tuning price. For the same experimental constraints 
(in particular the same lower limit on Mh), the minimal value of Ao is for tan/3 = 1.65 a 
factor of ~ 3 — 4 smaller than in JQ] (remember the factor 2 in the present definition of Ao), 
in agreement with previous findings ®0- For tan/3 ~ 0(10), one-loop corrections give 
much smaller effects, with typically a 30% reduction. In Figs. 1 to 3 we observe a very strong 
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Figure 1: 27ie price of fine tuning for tan (3 = 1.65 , as a function of various variables in the 
minimal supergravity model. An upper limit of 1.2 TeV on the heavier stop mass is imposed 
in the scanning. All experimental constraints described in the text are included. In all plots, 
except for A 0 versus M h , the bound M h >90 GeV is included. The mass of the lighter physical 
chargino and of the heavier physical stop are denoted by mci and M^ 2 , respectively. 
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Figure 2: As in Fig. 1, but for tanp = 2.5 




















Figure 3: As in Fig. 1, but for tan j3 
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dependence of A 0 on M h} which has been explained in Section 2. In consequence, the new 
lower limit M h > 90 GeV |12j pushes, for tan/? = 1.65, the minimal value of A 0 into the range 
~ 0(100) (~ 0(200) with the definition of A 0 used in ||). We note the increase by a factor 3 
in the minimal Ao with the change in the lower bound on Mh from 80 to 90 GeV. The results 
in all three Figures are qualitatively similar except for the overall decrease in the fine-tuning 
price with increasing tan /?. One more difference is that the n < 0 branch of solutions has 
disappeared at tan/? = 1.65. For so small a value of tan/?, as explained earlier, negative ji is 
no longer compatible with the present bound Mh > 90 GeV. 


The fine tuning decreases with increasing tan/?, with values of Ao marginally reaching 
Ao ~ 10 for 3 < tan/? < 15. This is shown in Fig. 4a, where we plot (solid line) the minimal 
Ao as a function of tan/? for Mh > 90 GeV. Also in Fig. 4a we show similar plots, but for 
hypothetical lower limits on the Higgs boson mass Mh > 100, 105, 110 and 115 GeV. We notice 
that for Mh > 115 GeV the fine-tuning is large for all values of tan/? < 30. It is also interesting 
to observe that the bulk of the parameter range shown in Figs. 1 to 3 gives interestingly large 
fine tuning, even for intermediate values of tan/?. Finally, given the striking dependence of A 0 
on Mh, it is interesting to see its dependence on tan/? under the assumption that we know the 
Higgs boson mass. In Fig. 4b this dependence is plotted for the hypothetical values M h = 95, 
100, 105, 110 and 115 GeV. For values Mh < 105 GeV, A 0 as a function of tan/? has a clear 
minimum, which moves towards larger values of tan/? with increasing Mh- 


4 Bottom-Tau Yukawa Unification and b —>■ sy Decay 


Up to now, we have been discussing fine tuning in the minimal supergravity model, without 
any constraints imposed on Yukawa couplings at the GUT scale. One important remark is that, 
in such a framework, the b —> sq decay, although constraining for the parameter space, does 
not have any impact on the necessary amount of fine tuning. The results for the minimal A 0 
presented in Figs. 1 to 3 and 4a does not depend at all on the inclusion of b —► sq decay among 
our experimental constraints for tan/? < 15, and are negligibly modified for tan/? up to 30. 
This is no longer true if we impose some constraints on the Yukawa sector, which is discussed 
in this Section and Section 7. 


One interesting possibility is b — t Yukawa-coupling unification at the GUT scale 14]. It is 


well known that exact b—r Yukawa-coupling unification, at the level of two-loop renormalization 
group equations for the running from the GUT scale down to Mz, supplemented by three-loop 
QCD running down to the scale Mb of the pole mass and finite two-loop QCD corrections at 
this scale, is possible only for very small or very large values of tan /?. This is due to the fact 
that renormalization of the 5-quark mass by strong interactions is too strong, and has to be 
partly compensated by a large t-quark Yukawa coupling. This result is shown in Fig. 5a. We 
compare there the running mass rrib(Mz) obtained by the running down from Mqut , where we 
take Yb = Y T , with the range of rrib(Mz) obtained from the pole mass Mb = (4.8 ±0.2) GeV ||34|| , 
taking into account the above-mentioned low-energy corrections. These translate the range of 
the pole mass: 4.6 < Mb < 5.0 GeV into the following range of the running mass rrib(Mz)- 
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Figure 4: Fine-tuning measures as functions of tan/3. In panels (a),(c) and (e), lower limits on 
the Higgs boson mass of 90 GeV (solid), 100 GeV (long-dashed), 105 GeV (dashed) 110 GeV 
(dotted) and 115 GeV (dot-dashed) have been assumed. In panels (b), (d) and (f), M h has been 
fixed to 95 GeV (solid), 100 GeV (long-dashed), 105 GeV (dashed) 110 GeV (dotted) and 115 
GeV (dot-dashed). Panels (a) and (b) correspond to independent Mi/ 2 , A 0 and ft parameters. 
In panels (c), (d) and (e), (f), linear dependences Mi/ 2 = o and A 0 = ca^F o, respectively, 
have been assumed. 
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Figure 5: a) The running mass rn b (M z ) obtained from strict b — r Yukawa coupling unification 
at Mqut = 2 x 10 16 GeV for different values of a s (M z ), before inclusion of one-loop super- 
symmetric corrections, b) The minimal departure from Y b = Y T at Mqut measured by the ratio 
Y b /Y T — 1, which is necessary for obtaining the correct b mass in the minimal supergravity model 
with one-loop supersymmetric corrections included. 


2.72 < m b (M z ) < 3.16 GeV. To remain conservative, we use a s (M z ) = 0.115(0.121) to obtain 
an upper (lower) limit on m b (M z ). 

It is also well known Ji5[ ^6| that, at least for large values of tan/?, supersymmetric finite 
one-loop corrections (neglected in Fig. 5a) are very important. These corrections are usually 
not considered for intermediate values of tan/? but, as we shall demonstrate, they are also very 
important there and make b — r unification viable in much larger range of tan /? than generally 
believed (see also ]37|). However, one has then to pay a higher fine-tuning price! 


One-loop diagrams with bottom squark-gluino and top squark-chargino loops make a con¬ 
tribution to the bottom-quark mass which is proportional to tan/3 ||35, | 


We recall that, 


to a good approximation, the one-loop correction to the bottom quark mass is given by the 
expression: 
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Figure 6 : One-loop supersymmetric corrections to the b-quark mass as functions of various 
parameters for tan (3 = 10 (left panels) and 30 (right panels), assuming the minimal supergravity 
scenario and imposing all experimental cuts except for the b —» S 7 constraint. Acceptable values 
of Mb are obtained for Am^/nib < —0.14 for tan (J — 10 and A mb/mb < —0.1 for tan (3 — 30, 
corresponding to the regions below the solid horizontal lines. 
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and the function I(a,b,c ) is always positive and approximately inversely proportional to its 
largest argument. This is the correction to the running rn b (M z ). It is clear from Fig. 5a that 
for b — t unification in the intermediate tan /3 region we need a negative correction of order 
(15-20)% for 3 ^ tan/3 ^ 20, and about a 10% correction for tan/3 = 30. According to ([TOD , 
such corrections require fi < 0, and their dependence on some other parameters is shown in 
Fig. 6. 

We notice that, as expected from m, b — t unification is easier for tan /3 = 30 than for 
tan /3 ~ 10. In the latter case it requires A t ^ 0, in order to obtain an enhancement in (|T0D or at 


least to avoid any cancellation between the two terms in ([id]) . This is a strong constraint on the 
parameter space. Since A t is given by (|9]), b—r unification requires large positive A 0 and not too 
large a Mg (i.e., M\/ 2 )- In addition, the low-energy value of A t is then always relatively small, 
and this explains the stronger upper bound on Mh seen in Fig. 6 (for a similar conclusion, see 

0 


We see in Fig. 5b that, for tan/3 < 10, the possibility of exact b — r unification evaporates 
quite quickly, with a non-unification window for 2 < tan < 8 — 10, depending on the value of 
a s . However, we also see that supersymmetric one-loop corrections are large enough to assure 
unification within 10% in almost the whole range of small and intermediate tan /3. 


For tan/3 > 10, the qualitative picture changes gradually. The overall factor of tan/3, on 
the one hand, and the need for smaller corrections, on the other hand, lead to the situation 
where a partial cancellation of the two terms in (]l^) is necessary, or both corrections must be 
suppressed by sufficiently heavy squark masses. Therefore, as seen in Fig. 6 , b — r unification 
for tan/3 = 30 typically requires a negative value of A t , and is only marginally possible for 
positive A t , for heavy enough squarks. A similar but more extreme situation occurs for very 
large tan /3 values, which will be discussed in Section 7. It is worth recalling already here that 
the second term in (|T0|) is typically at most of order of (20-30)% of the first term pjfl , due to ([)!). 
Thus, cancellation of the two terms is limited, and for very large tan /3 the contribution of (|T0| ) 
must be anyway suppressed by requiring heavy squarks. This trend is visible in Fig. 7a already 
for tan/3 = 30. The Higgs-boson mass is not constrained by b — t unification, since A t can be 
negative and large. Finally, we observe that, for any value of tan/3, the one-loop correction 
to m b (M z ) (|ig) remains approximately constant after simultaneous rescaling of /i, M%/ 2 and 
mo- Since proper electroweak breaking correlates /i with M\n and mo, the loop correction to 
m b (M z ) is weakly dependent on sparticle masses. 


Returning now to the fine-tuning price, we show in Fig. 7b the dependence of Ao on tan/3, 
with exact b — r unification imposed as an additional constraint on the parameter space. This 
dependence is shown both without and with b —> sy decay included among the experimental 
constraints, and we first focus on the case with b —> sq excluded. A comparison with Fig. 3 
(where b — r unification was not imposed) shows a substantial increase in the fine-tuning price 
for tan/3 = 10. This follows from the large values of Ao needed in this case for 5 — r unification 
(see the strong dependence of A 0 on A 0 in Fig. 3) and from the simultaneous ease in satisfying 
the b —> sy constraint for tan/3 ~ 10 (as will be discussed shortly). On the other hand, for 
tan/3 = 30 it is easy to have b — r unification. Hence, the price A 0 to be paid without imposing 
the b —> sy constraint in Fig. 7b (dashed line) is essentially the same with or without b — r 
unification (actually, it is slightly below the value seen in Fig. 4a for the case without b — r 
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Figure 7: a) Lower limits on the lighter (dotted lines) and heavier (solid lines) stop and on 
the CP—odd Higgs boson A 0 (dashed lines) in the minimal supergravity scenario with b — r 
Yukawa coupling unification, as functions o/tan/3. Upper (lower) lines refer to the case with 
the b —>■ sy constraint imposed (not imposed), b) The corresponding fine-tuning price with the 
b —> s'y constraint imposed (solid line) and not imposed (dashed line). 


unification but with the b —» sy constraint included). 


We turn our attention now to a deeper understanding of the b —> sy constraint and its 
interplay with b — r unification. The first point we would like to make is that b —> sy decay 
is a rigid constraint in the minimal supergravity model, but is only an optional one for the 
general low-energy effective MSSM. Its inclusion depends on the strong assumption that the 
stop-chargino-strange quark mixing angle is the same as the CKM element Vt s . This is the case 
only if squark mass matrices are diagonal in the super-KM basis, which is realized, for instance, 
in the minimal supergravity model. However, for the right-handed up-squark sector such an 
assumption is not imposed upon us by FCNC processes |38|]. Indeed, aligning the squark flavour 


basis with that of the quarks, the up-type squark right-handed flavour off-diagonal mass squared 
matrix elements and (niU)f( ri are unconstrained by other FCNC processes. Therefore, 

in the limit that the other flavour off-diagonal matrix elements are zero, and for sufficiently 
small the couplings of charginos to stops and strange quark read 
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= —hhVi'Z 2 ? sin#,- 


where Zj? are matrices diagonalizing the chargino mass matrix (defined in |32|), the h t ^, c are 


Yukawa couplings and the U are stop mass eigenstates: t R = — sin 6^2 + cos 6 ^ 1 . The factor 
( A R m 2 )23 can be considered as a free parameter of the low-energy MSSM. Indeed, there exist 
GUT models |I3| that predict the mixing factor in the vertex stkC'f to be considerably different 
from the CKM matrix element Vt s ■ We conclude that only a small departure from the minimal 
supergravity model is sufficient to relax the b —> sy constraint, and it is interesting to study 
separately its impact on the fine-tuning price. 

In the minimal supergravity model the dominant contributions to b —> sy decay come from 
the chargino-stop and charged Higgs-boson/top-quark loops. For intermediate and large tan/3, 
one can estimate these using the formulae of |28| in the approximation of no mixing between 
the gaugino and higgsinos, i.e., for M w -C max(M 2l |//|). We get 
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Ah+ 


A, 
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/| 7 3 m t 1 ) 

‘'Ml , r r) I 



(13) 

(14) 


(15) 


where ti(^ 2 ) denotes the lighter (heavier) stop, 

cos 2 e ; = 1 (1 + , a = - , A} = G F ~J a /(2Tf ViV a (16) 

h ti 


and the functions f^ k \x) given in |28[ are negative. The contribution Ac is effectively propor¬ 
tional to the stop mixing parameter A tl and the sign of Ac relative to Aw and Ah+ is negative 
for A t fji < 0 . 


We can discuss now the interplay of the b — r unification and b —> sy constraints. The 
chargino-loop contribution (|I5| ) has to be small or positive, since the Standard Model contri¬ 
bution and the charged Higgs-boson exchange (both negative) leave little room for additional 
constructive contributions. Hence, one generically needs A t [i < 0. Since fi < 0 for b — r unifi¬ 
cation, both constraints together require A t > 0. This is in line with our earlier results for the 
proper correction to the b mass for tan (3 < 10, [] but typically in conflict with such corrections 
for larger values of tan /3. In the latter case, both constraints can be satihed only at the expense 

lr This does not constrain the parameter space more than b — r unification itself. Note also that, if we do not 
insist on b — r unification, the b —> sy constraint is easily satisfied since /i > 0 is possible. 
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of heavy squarks (to suppress a positive A t correction to the 6 -quark mass or a negative A t 
correction to b —> sq) and a heavy pseudoscalar A 0 . Hence we have to pay a higher fine-tuning 
price, as seen in Figs. 7a,b. 


5 Linear Relations between MSSM Parameters 


The minimal supergravity model with universal soft mass parameters discussed so far is based 
on the assumption that scalar mass parameters are not independent of each other (and similarly 
for gaugino masses). This has obvious implications for the question of fine tuning, which can 
only be considered once the set of initial parameters is specified. In particular, one could relax 
the universality assumption and study the question of fine tuning for each sfcrmion flavour 
separately, with A^. being a measure of the fine tuning.. An increase (decrease) in the number 
of initial parameters is not directly correlated with increase or decrease of the necessary fine 
tuning. For instance, suppose the parameters aj and aj with derivatives A a . and A a . are 
assumed to be not independent but linearly related: a* = Cijdj. In this new scenario, the 
fine tuning is measured by A aiaj = A ai + A aj . The relative magnitudes of A aiaj , A ai and A aj 
depend on the relative signs of A a . and A a ., and vary from one region of parameters to another. 
However, as observed in || and indicated in Section 2, the scalar sector has little impact on the 
overall fine tuning, since the derivatives A ^ are generically smaller than the other derivative^]. 
Therefore, scenarios with correlated or uncorrelated scalar masses have similar fine tuning. 

The discussion in Section 2 shows that most often the largest derivatives are A /10 , A Ml/2 and 
A^q and, moreover, in the phenomenologically relevant parameter space they are of opposite 
signs. Let us take as an example again small tan/?. We see in (|7|) that A w < 0, whereas 
Ami / 2 > 0 for ji > 0 , which is necessary to maximize Mh. Also, sign(AA 0 ) = —sign(/rA 0 ), 
so that for /i > 0, sign(A At 0 A^ 0 )=sign(A 0 ) and for negative A 0 the derivatives A M and A ^ 0 
are of opposite signs. Since negative Aq is necessary for maximizing Mh , one expects that, by 
assuming there is some theoretical reason why some of these parameters are not independent, 
one may significantly reduce the fine-tuning price. It is easy to study the simplest case of linear 
relations between these parameters. Linear relations are also enough to obtain substantial 
reductions in the fine tuning, as we find in our full one-loop numerical calculations. 

We obtain the biggest reduction in the fine-tuning price by treating My 2 and fi 0 , or A 0 and 
/i 0 , as linearly related to each other, and the best choice depends on the value of tan/?. This 
is shown in Fig. 8, where we plot A 0 versus A M ^ and A A ofM) for several values of tan/?. For 
low tan /? (close to the infrared quasi-fixed point) the best effect is obtained for the A 0 — 
correlation, with the minimal fine tuning decreasing by a factor ~ 3.5. In Figs. 4c and 4e we 
plot minimal values of A Ml/2 ^ and A_a 0#1 as functions of tan /? for several assumed limits on M h : 
M h > 90, 100, 105, 110 and 115 GeV, and in Figs. 4d and 4f we show similar plots but for 
M h = 95, 100, 105, 110 and 115 GeV. The strong price reductions are evident. We see that 

2 The important implications of relaxing universality for the Higgs-boson mass parameters in the large-tan/3 
scenario (see next Section) has a different origin: it helps to permit electroweak symmetry breaking in a larger 
part of the parameter space. 
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Figure 8: Fine-tuning for correlated GUT-scale parameters (Mi/ 2 = cm^P o m the left panels 
and A 0 = ca^Po in the right ones) versus fine tuning for uncorrelated GUT-scale parameters 
for several values of tan /3. Universal soft scalar masses are assumed at the GUT scale. 
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Figure 9: As in Fig. 1, but assuming the linear correlation A 0 = o at the GUT scale 
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Figure 12: Correlations between the GUT scale parameters for several values of tan /3, assuming 
the following cuts on fine tuning. For tan/3 = 1.65: in the left panel, parameter sets with 
A Mu < 50 are indicated by stars, and sets with A M/i < 100 by points, in the right panel, 
sets with A AfJ < 20 are indicated by by stars, A^ < 30 by circles and A^ < 50 by points. 
For tan/3 = 2.5 and 5: sets with A M/X , A^ < 10 are indicated by stars and A M/X , A^ M < 30 
by points. For tan/3 = 30: in the left panel, sets with A Mfl < 3 are indicated by stars and 
A M/l < 10 by points, and in the right panel sets with Aa^ < 30 are indicated by points. 
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Figure 13: Upper limit on the lighter stop (a), stau (b), chargino (c) and neutralino (d) as 
functions of tan fj for universal soft scalar masses at the GUT scale requiring the fine-tuning 
measures A 0 (solid), A Mm (dashed) and A^ M (dotted) smaller than 10. 
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A ~ 0(1) is compatible with a large range of tan/3 values and Higgs boson masses. In Figs. 9, 
10 and 11, we plot A Ml/2lm and A a 0/j , 0 as functions of various mass parameters or physical 
masses, for several values of tan/3. The overall pattern remains similar to the A 0 case, but 
with an order of magnitude or more rescaling in the absolute values of the A’s. Nevertheless, 
putting some upper bound on the acceptable fine tuning remains a strong constraint on the 
superparticlc spectrum. This is more clearly seen in Fig. 13, where we plot the upper bounds 
on several physical masses as a function of tan/3, as obtained by requiring A 0 < 10. In this plot 
we compare the bounds obtained for M\/ 2 — po an d A 0 — p 0 correlations with the bounds for the 
uncorrelated case (A 0 < 10). As is seen clearly in Fig. 13, the upper bounds are considerably 
relaxed if correlations are imposed, suggesting that it is premature to use the fine-tuning price to 
derive convincingly any upper mass limits, in the absence of deeper theoretical understanding. 


6 String-Inspired Models 


So far, we have discussed linear dependences among soft supersymmetry-breaking terms in a 
model-independent way. We now take a more theoretical viewpoint, according to which the 
soft supersymmetry breaking parameters are predicted by some physics at the GUT or string 
scale. It is likely that they will emerge from the high-scale theory described in terms of more 
fundamental parameters. It is also plausible that the number of these parameters, at least 
of the relevant ones, is smaller than the number of soft terms. The latter will then not be 
independent. Such scenarios indeed emerge in various toy super gravity/string models for soft 
terms. The fine-tuning criterion would then require some revision: even if the number of new 
parameters is not smaller than the number of soft terms discussed earlier, a reparametrization 
may introduce more “natural” fundamental parameters. Generically, in supergravity models, 
one can write the soft terms as 

ch = m 3/2 fi (p a ) (17) 

where a, is one of the soft supersymmetry-breaking terms (M 1 / 2 , m 0 , A 0 , B 0 ), or /i 0 , and 
the gravitino mass m 3 / 2 sets the overall mass scale. The functions fi(p a ) are functions of 
dimensionless parameters p a which can be regarded as, e.g., angles determining the goldstino 
direction in the dilaton and moduli held space. 

Among the questions one can ask are: 

a) Does there exist a set of parameters p a which is more natural than the soft terms them¬ 
selves, and what are their properties? 

b) Are there any simple models with such parametrizations? 

Within such a framework, we should study fine tuning with respect to the parameters rn 3 / 2 and 
p a . In the case of m 3 / 2 , simple dimensional analysis tells us that 

A m , /a = 1 (18) 
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On the other hand, the general formula for p a is: 
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As an example, we now study a simple toy model [[40| in which soft supersymmetry-breaking 
terms and the Higgs mixing parameter /j,q are described at the GUT scale by the following 
parametrization: 


M 1/2 = V / 3m 3 /2 sin 9e l ' ys 
A 0 = —V3m 3 /2 sin 9e lls 

KJo = KJo = m ^ 2 ( X _ 3 cos2 9 ( 0 3 + 0 e)) 

Po = m 3 /2 (l + V3cos9 (o 3 e* 73 + 0 6 e* 76 )) 

B 0 p 0 = 2m 3 / 2 (l + \/3 cos 9 (0 3 cos y 3 + 0 6 cos 7 6 ) + 3 cos 2 9 cos (y 3 - 7 6 ) 0 3 0 6 


( 20 ) 


where the angles 9 , 0/ determine the goldstino direction in the dilaton/moduli parameter space, 
with sin# « 1 (0) corresponding to dilaton- (moduli-) dominated supersymmetry breaking, 
and the 7/ are phases which we set to zero for simplicity in the rest of this discussion. 


We now re-examine fine tuning in this new parametrization, considering first the sensitivity 
of Mf to sin0. We obtain A sin e from (|T9D using the soft parameters (|20| ) and the coefficients 
c'l obtained by solving the one-loop renormalization-group equations. The resulting formula is 
very complicated (remember that Mz, Ma and tan/3 in (pT|) depend on the parameters p a ) and 
will not be given here. It is, however, not very difficult to check that A S i n e has typically several 
zeroes as a function of 9 (for example, in the limit y —>■ 1 there are zeros for 9 = nir/2). Thus, 
there are regions in the new parameter space where the sensitivity to sind is small. However, 
this does not mean yet that one can easily avoid fine tuning, since it is necessary to check 
whether the regions of small A sin e correspond to phenomenologically acceptable solutions. To 
do this, we analyse Mf itself as a function of 9. Fig. 14 shows some typical plots of M^{9). 


We see that M^(9) is either negative or rather large (in units of 777 2 / 2 ) and positive at the 
extrema, and hence not acceptable. Experimental lower bounds on the masses of superpartners 
imply that the scale of supersymmetry breaking measured by m 3 /2 must be rather big compared 
to the weak scale. From a phenomenological point of view, the interesting regions of the pa¬ 
rameter space are only those which give positive but rather small values of M%. Unfortunately, 
A sin 0 is never very small in such regions. 


We have checked this by a numerical calculation. We have scanned the (0, 0 3 , @e) parameter 
space looking for solutions with Mf between 0 and and with Ma > 0.6 M Z - Such solutions 
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Figure 14: Dependence of M z on the angle 9 in the string parametrization of the soft terms 
shown in (20) for @ 3 = @ 6 = —0.5 (solid lines) and @ 3 = —0 6 = —0.5 (dashed lines). 

exist only in quite a small part of the (0 3 , 0 6 ) parameter space. Moreover, they give very small 
values of tan/3 (quite close to 1), and have A sin e always well above 100. Thus, we conclude 
that the parametrization (12^) cannot solve the fine-tuning problem. 

A search for more attractive models for soft terms, perhaps guided by the phenomenological 
discussion of Section 5, is certainly very important. 


7 Large-tan [3 Region 


In this Section we update the status of scenarios with (at least approximate) t — b — r Yukawa 
coupling unification |44], [42| and discuss their fine-tuning aspects. Such a possibility is realized, 
for instance, in S'O(IO)— type models. For m t = 175 GeV, t — b — r unification predicts large 
values of tan/3: tan/3 ~ 50, and the Higgs boson mass Mh ^110 GeV. Clearly, if the Higgs 
boson is not found at LEP 2, the phenomenological relevance of the large tan /3 region will be 
accentuated. 


Phenomenological properties of the large tan /3 region are well understood |T0 


Important aspects are the breaking of the electroweak symmetry and supersymmetric one-loop 
corrections to the bottom quark mass and to b —> sy decay. To organize our discussion, let 
us begin with exact t — b — r unification of Yukawa couplings in the minimal supergravity 
model. The results of Section 5 can be readily used to conclude that it is not a realistic 
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scenario 


36]. It is sufficient to observe that, in the parameter space constrained by requiring 


proper electroweak symmetry breaking, it is impossible to obtain sufficiently small one-loop 
supersymmetric corrections to the 6-quark mass for stop masses up to (9(10 TeV) (as can be 
estimated from (p~0|) . The problem is even worse if we try to be consistent with 6 —» sy decay. 


The question of some interest is how far we have to depart from exact unification of all 
three couplings in the minimal supergravity model to obtain a more realistic parameter space. 
One way of answering this question is to impose 6 — r unification, and to study the minimal 
values of the stop masses and of the fine-tuning measure Ao which are necessary to satisfy all 
the remaining constraints, as a function of tan/3. This is shown in Fig. 7a and 7b, respectively, 
requiring the correct value of the 6-quark mass and, optionally, the correct BR(b —> sy). As 
we discussed in Section 4, the prediction for the latter can be modified by a departure from 
the minimal supergravity model that admits some flavour structure in the stop mass matrices. 
We see from Fig. 7a that correct M& requires M f ~ 2 ~ 0.75(1.7) TeV for tan/3=45 without (with) 
the 6 —> sy constraint included. The corresponding values of Ao 0 are 170 and 35 respectively. 
For m t = 175 GeV, tan/3=45 corresponds to Y t /Yb ~ 1.6 (for My ~ 1 TeV). The origin of all 
these results is the extremely constraining role played by electroweak symmetry breaking in 
the minimal supergravity model with Y t ss Y h . In the limit Y t = Yb, the two soft Higgs boson 
masses rn 2 H} and m 2 H2 run almost in parallel, and electroweak symmetry breaking occurs only 
for very large values of My 2 and p. 


It has been pointed out in |44| that a qualitatively new situation appears in the large tan/3 

This 


scenario if we relax the universality of the Higgs-doublet soft mass parameters 45, 46, 44 


is because the correlation p M]/ 2 is no longer necessary for proper electroweak symmetry 
breaking, and one obtains solutions with p ~ M \/ 2 ~ (9(M^): as discussed in [Q, the hierachy 
2 H2 > rriQ is necessary for this. In consequence, in this scenario the supersymmetric 


m 2 H i m 2 


loop corrections ( |TC)| ) to irib(Mz) can be small. It is, therefore, interesting to repeat the analysis 
in this case. Since it is easy to obtain acceptable physical Mb , even for Y t = Yb = Y T , we restrict 
our analysis to this case. We study the case tan/3 = 50 and impose Mb = 4.8 ± 0.2 GeV, i.e., 
2.72 — 3.16 GeV for mb(Mz). 


In Fig. 15 we show some results for A 0 as a function of several mass parameters, with all 
the constraints included except for 6 —> sy. Only p < 0 is possible since, as is clear from 
Fig. 5a, only negative one-loop supersymmetric corrections are compatible with the correct 
bottom-quark mass. Moreover, the correction has to be small enough and, therefore, A t tends 
to be negative and squarks must be relatively heavy. Due to the hierachy m 2 H /$> tUq, the 
pseudoscalar A 0 is heavy enough to assure a small amount of foie tuning A 0 ~ 10 is possible. 
This result makes the large tan /3 region quite acceptable from the naturalness point of view. 


If we insist on being consistent with 6 —> sy decay, the fine-tuning price increases to A > 40, 
as seen in Fig. 16. This happens for reasons similar to those discussed for tan/3 = 30 in the 
minimal model, and the discussion at the end of Section 4 applies unchanged to the present 

3 For large tan/3 it is important to consider the derivatives of Mz and tan/3, since the latter are proportional 
to tan/3 and can be large. 

4 Typically the dominant derivatives are ( di/ tan /3)(d tan /3/dcii) « — tan /3((Bc^ai + ^c l B ai)/M\). 
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Figure 15: The price of fine tuning for tan (3 = 50 and t — b — r Yukawa coupling unification, as 
a function of various variables in models with non-universal Higgs boson masses. The b sj 
constraint is not included. 
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case. 


Finally, we remark that the non-universal Higgs boson masses discussed here give solutions 
with higgsino-like neutralinos. Thus, such scenarios generically lead to a rather low neutralino 
dark matter density |3(| . 


8 Conclusions 


Comparing the situation before and after LEP, the fine-tuning price in the minimal supergravity 
model has increased significantly, largely as a result of the unsuccessful Higgs boson search. 
Comparing different values of tan/3, we find that naturalness favours an intermediate range. 
Fine tuning increases for small values because of the lower limit on the Higgs mass, in particular, 
and increases for large values because of the difficulty in assuring correct electroweak symmetry 
breaking. 

Additional theoretical assumptions may have a significant impact on the fine-tuning price. 
For example, requiring b—r Yukawa-coupling unification would increase the price significantly at 
intermediate tan/3, whereas imposing certain linear correlations between mass parameters could 
diminish it substantially. One particular class of models imposing such have been motivated 
from string constructions: unfortunately, those currently available do not seem to reduce the 
fine-tuning price significantly, so naturalness considerations do not favour these models to any 
substantial extent. However, the search for realistic theoretical models which do reduce the 
fine-tuning price is a very interesting issue. 

We have found that b —> sy decay is a potentially important constraint for large tan/3, 
but we would argue that it should be regarded as optional. The flavour structure of squark 
couplings could differ from those of the quarks, and there are no direct FCNC limits on flavour 
violation among superpartners of the up quarks. 

A final comment concerns the region of very large tan /3. In this case, the fine-tuning price 
can be reduced quite substantially by allowing non-universal soft mass parameters for the Higgs 
bosons. This is in contrast to the situation at lower tan/3, where non-universal mass parameters 
do not reduce the price significantly. 

We re-emphasize that naturalness is subjective criterion, based on physical intuition rather 
than mathematical rigour. Nevertheless, it may serve as an important guideline that offers 
some discrimination between different theoretical models and assumptions. As such, it may 
indicate which domains of parameter space are to be preferred. However, one should be very 
careful in using it to set any absolute upper bounds on the spectrum. We think it safer to use 
relative naturalness to compare different scenarios, as we have done in this paper. 
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